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Abstract. We study the Euclidean effective action per unit area and the charge 
density for a Dirac field in a two— dimensional spatial region, in the presence of 
a uniform magnetic field perpendicular to the 2D-plane, at finite temperature 
and density. In the limit of zero temperature we reproduce, after performing 
an adequate Lorentz boost, the Hall conductivity measured for different kinds of 
graphene samples, depending upon the phase choice in the fermionic determinant. 
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Graphene is a bidimensional array of carbon atoms, packed in a honeycomb 
crystal structure. Actually each layer of a graphene sample can be viewed as either 
an individual plane extracted from graphite, or else as an array of unrolled carbon 
nanotubes. Stable mono-, bi- and multi-layer samples of such material have been, 
recently and independently, obtained by two groups [TJ [2j , and a surprising behavior 
of the Hall conductivity and related density of states in mono-layer samples has been 
unravelled. Even more recently, an equally unexpected, though different, behavior 
was reported [3] for bi-layer samples. 

The most remarkable feature of graphene's structure, from the theoretical point 
of view, is that its quanta, or quasi-particles behave as two species (to account for 
the spin of the elementary non-relativistic constituents) of massless relativistic Dirac 
particles in the two non-equivalent representations of the Clifford algebra (which 
correspond to the two non-equivalent vertices in the first Brillouin zone [HH]) with 
an effective "speed of light" about two orders of magnitude smaller than c . 

To the best of our knowledge, the first approach to the quantum Hall effect in 
terms of a quantum relativistic Dirac field theory, at finite temperature and chemical 
potential, appeared in reference [6], where a dimensional reduction argument was used. 
Later on, the relativistic Hall conductivity was obtained in 7, 8 using Green's function 
methods (for an entirely different calculation see, e.g., [9]). In references [6l[71[8], 
the divergent series were regularized through definitions which, as we will explain, 
are equivalent to neglecting the phase of the determinant. In the zero temperature 
limit, the results of refs. [3 OH] reproduce the unexpected behavior of both the Hall 
conductivity and the density of states, as measured in mono-layer graphene [HE]. 

In the couple of papers pUHH], two of the authors of the present letter developed a 
finite temperature field theory calculation based upon (— function regularization of the 
Dirac determinant, and obtained the partition function and the related Hall current 
and density of states. There, the phase of the determinant was included, and its sign 
fixed according to the conventional wisdom [12], which lead to a Hall conductivity 
displaying a plateau around zero chemical potential. 

It is the aim of this paper to show that, in turn, the inclusion of the phase of 
determinant with opposite sign leads to a Hall conductivity and to a zero-temperature 
density of states which coincide with the ones recently reported in the case of bi-layer 
graphene [3] , while the behavior of mono-layer graphene is reproduced when the phase 
is ignored, in coincidence with the results in [6l[7l[8j. 

In order to study the temperature-dependent effects for the system at hand, we 
will consider the three-dimensional (3D) Euclidean space, with metric (+,+,+) and 
coordinates {x,y,T = —it), the 7-matrices ji = o",, i = 1,2,3 and introduce the 
chemical potential as an imaginary component of the gauge potential [13] . Then, 
we will let the Euclidean imaginary time coordinate vary according to < r < /3, 
where (3 = l/ksT, ks is the Boltzmann constant, and impose antiperiodic boundary 
conditions on the Dirac field to reproduce Fermi-Dirac statistics. As it is well known, 
another faithful non-equivalent representation of the Clifford algebra exists in odd 
dimensions, in which one of the gamma matrices changes sign (or, equivalently, all 
of them do so). We will comment about the consequences of such a change of 
representation, wherever adequate, throughout the rest of the paper. 

Once some suitable regularization has been introduced, the partition function in 
the grand-canonical ensemble is given by 

In Z = In det(£>) = In dct(i/9+ eyl) , (1) 
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where — e is the electron charge. In order to evaluate the partition function in the 
zeta regularization approach [T3], we must determine the spectral resolution of the 
Euclidean Dirac operator, in the presence of a gauge potential = (0,Bx,i/j,/e) , 
which corresponds to the selection of a non-symmetric gauge for the magnetic field 
orthogonal to the plane - in what follows B > . Here below we just report the main 
results - for a detailed calculation see, e.g., [101 E]. The equation to be solved is 
(from here on natural units are used, i.e., H = c = 1, unless explicitly stated) 



After writing 



where 



and 



[<j\ id x + o~ 2 (id v + eBx) + cr 3 (id T + ifx) — u] \& = . 
Vk,l(x,y,r) = (2tt/3)~ 1/2 exp{iky + irXi}i(Jk,i(x) , 



<•,..,(..•) = j ^ M Si i - /- : 

1 Xk,i{x) ' 



Xi = (21 + 1)(tt/P), I el 



in order to satisfy the antiperiodic boundary conditions, one finds two types of 
eigenvalues and corresponding sets of eigenfunctions. 
[ I ] Asymmetric part of the spectrum 

loi = — (21 + l)n/[3 + ifi, l = — oo, ...,oo 

The corresponding eigenfunctions are quite particular, in that they are 
eigenfunctions of as with eigenvalue +1 , the ones in the orthogonal subspace being 
indeed eliminated by the square-integrability condition (for some related references 
see, for instance, [15 ). As a consequence, the corresponding eigenvalues uji are not 
the square roots of the eigenvalues of D'D. They will eventually lead to a spectral 
asymmetry and, thereby, to a phase of the determinant, which will be studied in detail 
below. 

[ II ] Symmetric part of the spectrum 



ui l n = ±y Xf + 2neB, n = 1, oo, I = — oo, oo 

For both kinds of eigenfunctions, the degeneracy per unit area is given by the 
well known Landau factor = eB/2-K . 

It is worthwhile to remark that, had we chosen the other non-equivalent 
representation of the 7-matrices, the eigenvalues of type [I] would have changed 
their sign, and the corresponding eigenfunctions would have been eigenfunctions of 173 
with eigenvalue —1. However, as will be discussed below, this fact will not lead to any 
modification in our physical predictions as long as /1 is real. Thus, considering the 
contributions of both noncquivalcnt representations will amount to an overall factor 
of two. 

When parity is defined as, e.g., in ref. [TB], it is easy to check that, for a general 
Dirac operator, the effect on the spectrum is uo p — —lo. This symmetry is obviously 
respected by the symmetric part [II] of our spectrum, while it actually produces a 
change in the sign of the asymmetric portion [I]. When acting on the latter, it is 
equivalent to fi — > — /i and, thus, to charge conjugation (ip(x) — > 72 ^(a:); A„(x) — > 
— A*(x); uj — > u! C = ui*). So, parity is broken, already at the classical level if only one 
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representation of the gamma matrices is considered, due to the square-integrability 
condition. 

If complex values of (i are allowed, the whole spectrum has an interesting 
symmetry: it turns out to be invariant under /i/3 — > n/3 + 2nik , k G Z , which is 
nothing but the symmetry under "large" gauge transformations. The conflict between 
this last symmetry and parity invariance in different regularization schemes is well 
known [17] , and created some controversy in the past [El HH [19] . 

In this letter, we will concentrate on the case of a real chemical potential fi. A 
discussion of both symmetries in the case of a complex chemical potential will be 
reported elsewhere |20| . 

Starting from the above described spectrum, we shall evaluate, according to (QJ 
the Euclidean effective action per unit area (in the statistical mechanics terminology, 
the latter coincides with the grand-potential per unit area in units of ksT). From this 
effective action, the mean fermionic number per unit area, Np, and, thus, the charge 
density can be retrieved as follows 

-idSeff 



S eff = logZ = In det(-D) . 



dfi 



(2) 



Here, the symbol = stands for the definition through an adequate regularization. 
For the reasons we have just explained (viz., all the eigenvalues are paired), it turns 
out that the contribution to the effective action coming from the symmetric part [ II ] 
of the spectrum does not suffer from regularization ambiguities. For instance, after a 
proper definition in terms of the ^-regularization 10J, it is given by 



oil 



C«(-l/2) (2eBfl 2 + £ log { (1 + * e"^") (l + z~ x e~^) } 

=1 



(3) 



with z :— exp{(3fj,} , e n = \j2neB . 

The contribution to the effective action arising from the asymmetric part [I] of 
the spectrum is given by the formal expression 



S eff 



1=0 



log{(-l)[(2/ + l)7r//3-z M ] x l(2l + l)ir/f3 + in}} 



(4) 



Choosing a symmetric regularization, as done in [U El [8] , is equivalent to ignoring the 
infinite term YaZq 1°s(~1) > which reduces the previous expression to 



Si ff = A L J2 log { [ (2/ + 1W/3] 2 + /i 2 } 



1=0 

i.e., one evaluates the logarithm of the "absolute value'' 
operator which, once regularized, leads to 



°eff 



A L log[2cosh(/x/3/2); 



of the Euclidean Dirac 



(5) 



However, a first-principle ^-function regularization of the determinant unavoidably 
drives to a careful definition of the phase of the determinant, which is equivalent to the 
selection of a cut in the complex plane of the eigenvalues [21] , when the asymmetric 
part of the spectrum is treated. Going back to the formal relation ([4]), we define it in 
proper mathematical sense in terms of the ^-function regularization: namely, 



E 

. (=0 



(2/ + l)-+z/i 



E 

1=0 



e ±m (2l+l)- + ifi 

P 
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the phase of the determinant being fixed by the cut in the complex plane of the 
eigenvalues [55]. More explicitly we can write 



E 

- 1=0 



(2/ + l)-+z/i 



E< 

1=0 



-is8 



(2l + l)-+ifj,e 

P 



-iO 



where —tt < 9 < 7r . The prescription usually adopted |12j amounts to choosing the 
cut in such a way that the expression in the last square bracket never vanishes, as 
one goes continuously from the eigenvalues with a positive real part to the eigenvalues 
with a negative real part, i.e., (21 + 1)tt//3 + /isin(9 and /icosB do not simultaneously 
vanish (which could happen if [i = (21 + l)ir//3 and 9 = —n/2 or fi = —(21 + l)n/(3 
and 9 = ir/2). This requires that the cut is chosen below (above) the real axis when \x 
is positive (negative), which gives for the final value 9 = 7rsign/i. With this choice, the 
contribution of the asymmetric part of the spectrum - see ref. |10j for more details - 
to the effective action is given by 



a eff 



1 13 log (2 cosh 



2 ' lr 1 V 2 
The opposite, and less popular, definition of the phase would lead to 



P\/j,\ +log 2 cosh 



M/3 



(G) 



(7) 



At this point, it is important to stress (always in the case of a real chemical 
potential) that exactly the same results are obtained, provided the same criteria are 
applied, if the other non-equivalent representation for the 7-matrices is chosen. Thus, 
the inclusion of both contributions, amounts to an overall factor of two, if the phase 
is consistently chosen in both representations. In this sense, the exclusion of the 
phase is equivalent to the adoption of opposite criteria for the phase selection in both 
representations. 

Putting together the contributions from the symmetric part (|3|) and the 
asymmetric part of the spectrum ([5]), J6J or ([7]), depending on the phase definition 
adopted, we come to the following expression for the (-function definition of the 
Euclidean effective action, 



5', 



eff 



A L < log 2 cosh 



HP 



. K /3H+/?G?(-1/2)V2^B 



log [(1 + ze~ Pen ) (1 + z 



n=l 



K = 0, ±1 , z — exp{/3 /i} , e n = \ 2neB . 
Here, k — corresponds to a vanishing phase, k = — 1 to the usual phase choice and 
K = +1 to the opposite and unusual phase choice. Note that e n is the absolute value 
of the n— th non-vanishing Landau level. In all cases, the Euclidean effective action 
is an even function of fi . Thus, it is invariant under charge conjugation and parity. 

Also in all the cases, the mean fermionic number per unit area, turns out to 
change sign under fi — > —fi . In fact, from its very definition in (|2|), we get 



1 , n/3 1 
-tanh — + -«;sgn(/z) 



71=1 



OO 

-E 

71=1 
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Note that the first two terms are the ones coming from the asymmetric part of the 
spectrum. So, no matter how one defines the phase of the determinant, at any finite 
temperature there is a kind of parity breaking charge, which is the sum of a /i-analytic 
contribution (the first term) and a non-analytic one (the phase of the determinant). 
In the zero temperature limit, for n < [i 2 /2eB < n + 1 , we finally obtain 

lim N{k;P,h) = (n+i^ ) sgn(M)Ai. (8) 
p^oo y z i 

In this limit, the contribution of the asymmetric part of the spectrum is the one 
corresponding to n — 0, and is non-analytic in all cases. It only vanishes, so that 
parity and charge conjugation symmetries are indeed fulfilled, for the most commonly 
accepted selection of the phase k = — 1 . 

Moreover, it is easy to check that Nernst's theorem holds true for any k . Actually, 
it turns out that the entropy can be obtained from the well known Boltzmann-von 
Neumann formula 

S((3, fi; B, k) = k B (l - P jj^j S e ff ; 

whence, one can verify by direct inspection that, indeed, 

Urn S(j3,fi;B,K) = 0, V/c = 0,±l, 

in agreement with Nernst's theorem. 

From ([5]), the mean value of the charge density in the zero temperature limit 
can be immediately obtained for one representation and one fermion species. Turning 
back to physical units, and recalling that the particle charge is — e , we find 

J t (k;/3,/x) = - eN F (K;/3,^) 

p^co , \ { 1 + k\ e 2 

— ► - sgn(^) 



he 

for n < (pi 2 /2eBhc 2 ) < n + 1 , n £ N , 

the spatial components of the current density being equal to zero in the absence of 
electric fields. Now, the zero temperature limit of the same vector in the presence of 
crossed homogeneous electric E' and magnetic B' fields can retrieved, for E' < B', by 
performing a Lorentz boost with absolute value of the velocity v = cE' / B 1 . Suppose, 
for definiteness, that the homogeneous electric field points toward the positive Oy- 
axis. Then, the speed of the Lorentz boost must point toward the negative Oy-axis 
and the transformation law gives as a result 

e 2 B' ( 1 + k\ . . 
J t {k;h) = \n-\ —J sgn(/z), 

e 2 E' ( 1 + k\ 

= — t— — 7T- w): 



for n < (/i 2 /2eBhc 2 ) < n + 1, n 6 N. As a consequence, the contribution to the 
quantized Hall conductivity at zero temperature becomes, for each representation and 
each fermion species, 

e 2 ( 1 



h 



sign(/i) , n £ 
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As explained throughout the paper, this result must be multiplied by an overall factor 
of four, in order to make contact with the relativistic effective theory associated to 
graphene [U[S]. It is interesting to remark that the three values of the phase correspond 
to three different vacuum polarization (Casimir-type) effects due to the interaction 
with the magnetic field at zero temperature. More precisely, the quantum (in the 
field-theory sense) filling factor is given, in each of the three cases, by 



The (rescaled) Hall conductivity is presented in figure [TJ for the three values of 
fij , as a function of vq — sgn(fi) fi 2 /2eBhc 2 , which is nothing but the classical density 
of carriers in the relativistic theory, divided by the total degeneracy of each Landau 
level. 

It is well known that any regularization procedure is acceptable, unless either it 
manifestly violates some of the symmetries of the system, or it is ruled out by the 
experimental data. So, the measured Hall conductivities of graphene as reported, for 
instance, in reference [3J, should shed light on the relevance (or lack thereof) of the 
different phases of the determinant. 

On these grounds, the first and clearer conclusion of this letter is that the behavior 
of monolayer graphene, as presented not only in [3] but also in [TJH], corresponds to 
k = 0, i.e., to not including the phase of the determinant, as done in [HI El [H] • In fact, 
in this case the (rescaled) Hall conductivity shows a jump of height 1 for vc — 0, and 
further jumps of the same magnitud for vc — ±1, ±2, .... 

Let us now compare our predictions with the contents of reference [3J, which is 
devoted to bilayer graphene. In this case, the (rescaled) Hall conductivity presents 
a jump of height 2 for vc — 0, and further jumps of height 1. The main point here 
concerns the positions of these subsequent jumps. As a matter of fact, according to 
figure l.b in the same reference, these subsequent jumps appear for vc = ±1,±2, 
which is exactly the behavior predicted, in our calculation, for k = + 1 (the less 
popular selection of phase in the Dirac determinant). However, the same reference 
interprets the Hall behavior of bilayer graphene through a theoretical prediction first 
made in [23j . where the theory is "almost" nonrelativistic, with a Landau spectrum 
given by E n = ±^^-y / n(n — 1). This last model does, indeed, predict a double jump 
for vc = 0, due to the existence of two zero modes in each representation. But, for the 
very same reason, the next jump should appear, in this theoretical scenario, at vc = 
±y/2 ~ ±| (for a related discussion see, for instance, [24]). As stressed by the authors, 
figure l.b in [3] is only schematic. However, the experimental results corresponding 
to B — 12T in figures 2.b and 2.c of the same reference also tend to confirm the 
prediction of our relativistic quantum field calculation, with k — 1 (unusual phase), 
where the next jumps occurs at vc = ±L The text in the same reference also seems 
to confirm our prediction, since the distance between jumps is said to be, always for 
B = 12T, An ~ 1.2 x 10 12 cm -2 , which corresponds to Ave = An-J^g ~ 1, the same 
for all jumps. However, the experimental results corresponding to B — 20T seem to 
agree with the width of the first plateau being approximately 1.3 times the width of 
the subsequent ones. Thus, further experimental study of bilayer graphene is crucial 
in distinguishing between both theoretical scenarios. 

At this point, one can naturally wonder whether there is place at all for the 
usual selection of phase (n = —1) in the description of graphene samples. It is quite 




(9) 
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Figure 1. Hall conductivity for different selections of the phase of the 
determinant. Top to bottom: k = 0,k = 1,k = —1 . In all cases, the horizontal 
axis represents up = sgn(/x) /2eBHc 2 and the vertical one, a xy h/4e 2 . 



interesting to gather that the three non-equivalent phase selections correspond to 
the three non-equivalent unitary representations of the cyclic group C3, which is 
precisely the relevant symmetry group for graphene. Even though the study of the 
Hall conductivity in graphene samples with three layers [25) is certainly far from being 
conclusive, they seem to indicate that a quantum Hall effect does occur in such devices, 
with a plateau at vq — 0. Does the usual phase selection correspond to the behavior 
of three-layered graphene? 

In any case, from a theoretical point of view, further experiments on graphene 
samples can give an answer to a long-standing question in the field of the zeta- 
function regularization, i.e.: Which phase must be selected in the definition of Dirac 
determinants, in order to evaluate effective actions? 
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